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Abstract 

The one-loop effective action for a slowly varying electromagnetic field is 
computed at finite temperature and density using a real-time formalism. 
We discuss the gauge invariance of the result. Corrections to the Debye 
mass from an electric field are computed at high temperature and high 
density. The effective coupling constant, defined from a purely electric 
weak-field expansion, behaves at high temperature very differently from 
the case of a magnetic field, and does not satisfy the renormalization group 
equation. The issue of pair production in the real-time formalism is dis- 
cussed and also its relevance for heavy-ion collisions. 
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1 Introduction 



The formulation of quantum field theory in an external field is interesting because of 
the many applications where the background field is strong |l|] and cannot be treated 
perturbatively. Some relativistic examples are the magnetic field in neutron stars, 
white dwarfs and heavy-ion collisions ||. More extreme situations are given by the 
electroweak phase transition and cosmic strings ||. Strong (colour) electric fields 
furthermore appear in some models of hadronization in heavy-ion collisions (for some 
recent discussions see e.g. [§J [H). In most of the examples above the fields exist in a 
thermal heat bath or some non-equilibrium background which is very different from the 
vacuum. In the case of hadronization in heavy-ion collisions it has actually been argued 
(see e.g. Ref. 0) that there is a time-interval during which local thermal equilibrium 
has been achieved but the external (colour) electric field has not yet been depleted due 
to particle production. For such a time-interval the results of the present paper apply. 

In this letter the one-loop effective action for a constant (or slowly varying) elec- 
tromagnetic field || is generalized to finite temperature and chemical potential, (T, /i). 
Or, to put it differently, the thermal effective action for a constant magnetic field f?j is 
generalized to arbitrary constant electromagnetic fields. The case of a pure magnetic 
field has been treated earlier in || |9], |7j and the incomplete result of || was corrected 
in [^] . The case of a general electromagnetic field has been studied in |llj and some 
corrections to the expressions for the effective action in those papers are presented here. 
For constant external fundamental fields there are various methods of calculating the 
effective action to all orders in the field. Schwinger calculated the effective action 
for a constant background field-strength in QED, for which the gauge fields are non- 
constant 0. The proper-time method used in || relies on the fact that the solution of 
the quantum mechanical equations of motion for a particle in the background field can 
be found explicitly, and an extension to finite (T, /x) is possible using real-time thermal 
propagators. The case of a magnetic field has a rather clear physical interpretation 
at finite (T, /x) with the particles in the heat bath occupying the time-independent 
Landau levels. This situation has been studied in detail in 0. In the presence of 
a slowly varying electric field the particles in equilibrium screen the field over a dis- 
tance determined by the Debye screening length. Schwinger's calculation for a general 
slowly varying electromagnetic field can be extended to finite (T, /i) in a rather straight- 
forward manner without finding the particle spectrum explicitly. The constant field 
approximation is satisfied if the relative gradient of the field is smaller than any other 
scale, i.e. \d a F^ v \/\F^ u \ <C \eF^ lu \ 1 ^ 2 , m e , T or n 1 / 3 , where n is the particle density. In 
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this letter we consider only the non-interacting e + e~-gas in a time-independent and 
slowly varying background where the conditions above are satisfied. For instance, the 
result applies to a shallow potential well with an arbitrary constant magnetic field and 
low density. It has been strongly argued against the possibility of having an electric 
field and a net charge density at equilibrium |12|]. We want to stress that this may be 
true for a strictly constant electric field, but it is perfectly physical to have a slowly 
varying field at thermal equilibrium and consider an expansion in derivatives in the 
field, at least for a stable system. A deeper analysis is required to determine whether 
the result can be used for more general situations. We shall limit the discussion here 
to the thermal corrections since the vacuum part is easily added when needed. 



2 One— loop effective action 

Schwinger's equation for the one-loop effective action, generalized to finite temperature 
and density, can be written as 

= *Tr(x| 1 - f F ( Po ,A) ( , * - , 1 . ) \x) , (1) 
am f — m + te \tf — m + it tf — m — it J 

where Tr is the trace over spin and x, if = 7 M (p jU — eA^), and fpipo, A ) is the thermal 
distribution function. Under the gauge transformation A^{x) \— > A^(x) + ^A(x), the 
zero temperature part of Eq. (|I|) is transformed to 

iTi(x\e- ieK{x) - e iehix) \x) , (2) 

and is thus gauge invariant. Since the thermal part of Eq.(|l]) is not invariant under a 
general gauge transformation we have to explain the apparent problem of gauge depen- 
dence. The density matrix for the whole system of fermions and the electromagnetic 
field is given by 

exp[-P f d*x(T™ D (x) -^ n 
PQED trex V [-Pld*x(T%%(x)-iie^n ' 1 J 

where Tq ED is the energy-momentum tensor 
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Tq ED = F» a F a v + (d a F» a )A u + im^d u m + -g l * u F a pF a P - g» u m(i$ - e4~ m)m . (4) 



One can easily verify that Tq ED is gauge invariant when Maxwell's equation d a F a ^ = 
e\&7 Al \I/ is satisfied. It is, however, common to consider the fermions in the background 
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field F^ v separately and only use the fermionic part of Pqed to define the density 
matrix. The energy defining p e + e -, 

P = Jd 3 xT^ = Jd 3 x$(iiidi - ejiAi + m)tf + e* 7o A ^] = j d 3 xUi lQ d<$> , (5) 

is then gauge dependent. It is not even a conserved quantity in a general gauge, in 
contrast to the total energy of the particles and the field. The equilibrium p e + e - can 
also be defined as the state with maximal entropy for a given expectation value of Po. 
It is then clear that the separation of the background field and the fermions is only 
meaningful in a gauge where Po is a conserved quantity. Only for time independent 
A^x) is Po separately conserved (assuming F^ u to be constant) and that determines 
the gauge we have to use, namely OqA^x) = 0, up to time independent gauge trans- 
formations which, anyway, leave the final result invariant .0 To be more precise, A(x) 
can a have time dependence of the form A(x) = c ■ t + A(x), where c is a constant 
and A(x) is an arbitrary function. This transformation shifts the potential Aq by the 
constant c. Such a constant is actually absorbed into the definition of the chemical 
potential, which is not an independent physical parameter but is to be determined by 
a given charge density. Only the difference p — eA is physically meaningful for p e + e -. 
Note that this discussion is equally relevant when = but then one normally puts 
A„ = Q. 

To find thermal one-particle distribution function fp{po, A ) we argue that at high 
temperature it must be reduced to the classical Boltzmann distribution for electrons 
and positrons, exp[— /?(Vp 2 + m ' 2 ± e A) -F A 4 )]) an d it is the sign of (p — eA ) that 
distinguishes between particles and anti-particles. We therefore write 

, / a \ 0(Po ~ e^o) , g(~ go + e^o) ,^ 

MPo, A)) - e/3(P0 _ M) - x + ef3{ . po+ll) + i ■ W 

The distribution function does not have to be chosen to represent an equilibrium dis- 
tribution. Other choices may be appropriate when the electric field drives the system 
out of equilibrium. It has been emphasized in the literature (see e.g. Ref.|p.2||) that the 
distribution in Eq.(^) has a non-trivial limit when T — > 0, and describes a Dirac sea 
filled up to the energy p. 

From Eq.(|H) one can follow the calculations in || to obtain the thermal part of 
the effective action. The only difference is that the trace should be taken in the basis 

3 Notice that such a gauge was not used in [10] . 
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\p , x )- We find 
dm 2 



where 
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2tt 3 / 2 

x(h(s) - ie)~ 1/2 exp 



30 dpo i f°° 

— — /f(po, v4 ) Im< / dse 2 ab cot (esa) coth(eso) 

00 27T JO 



—i[m 2 — ie)s + i 



(p - eA ) 2 



. 7T 
^ — 

h(s) -ie 4 



(7) 



h(s) 



2 j2 

a — 



ab 
B 



(eFcotheFs)oo , 
dfj,A u — d v A^ , 

B 2 — E 2 = IT 



- TP 

2 ^ 



E ■ B = ^ = -F^F**" 



IBI 



E = |E| 



(8) 



We have added — to h(s) in Eq.(0) in order to get the correct value using the principle 
branch when h(s) < 0. It is determined by the ^o-integration of Schwinger's formula. 
From the cot (esa) factor we find that the s-integral goes through a number of poles 
on the real axis, s = kir/ea, which were not apparent in the original expression. In the 
case of a pure magnetic field it has been shown that the poles are absent if we only 
include a finite number of Landau levels, and to get the correct result after summing 
over all Landau levels the s-integration contour has to go slightly below the real axis. 
We assume that the same prescription is valid even for non-zero E field since the poles 
are related to the B field. The expression in Eq.(]7|) can be directly integrated with 
respect to m 2 . The singularity at s = should be cured by subtracting the F^ v = 
part or using a ^-function regularization as in (or using Eq.flHSD). We notice that in 
the limit T — » 00, and fp(po, A Q ) — > 1/2, Eq.(|7]) equals the negative of the real part of 
the vacuum contribution up to divergent terms which are only quadratic in the field. 

It is possible to find an explicit expression for h(s) using the following observations. 
The characteristic equation for the fieldstrength tensor F = F u is 



det(AH - F) = A 4 + 2J^A 2 - Q 2 







with the eigenvalues 



\ = ±i(VT 2 + g 2 + F) 1/2 , ±(y/J r2 + g 2 -T) 



1/2 



(9) 



(10) 



We know that (F coth F) has a formal power series expansion involving only even pow- 
ers of F. The Cayley-Hamilton theorem states that F satisfies its own characteristic 
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equation, which can be used to reduce the powers of F down to at most F 2 , so the 
matrix structure must be 

FcothF = 7 11 + 5F 2 . (11) 

Taking the trace and determinant on both sides of Eq. ([TT|) gives two equations to 
determine 7 and 5. We find 



h(s) 

e 
2 



1/21 



+ I f 1 + 2^+g 2 ) ( v/ ^+^ ~ C ° th Mv^+g* - -^) 1/2 



(12) 



The relation s ■ /i(s) > is thus not valid for general electromagnetic field, contrary 
to the claim in [TTJ. The quantity ^=1,2,3 [ln(sinh(eFs)/eFs)] M , appearing in the 



effective action in [10], can be calculated using the same method, and we do not agree 
on their result either. The standard way to proceed from Eq.(|7]) is to deform the s— 
contour to the imaginary axis. Here, it is not always useful, except in special cases, 
since h(s) has zeros on the imaginary axis which give essential singularities from the 
exponential. On the other hand, the same problem occurs also on the real axis. To be 
more specific we can consider the following special cases 

I. E || B : h(s) = eEcothesE , 

B 2 



II. E1B, E > B : h(s) = — E 2 VE 2 - B 2 coth esVE 2 - B 2 , 

E 2 — B 2 \ es J 

III. E1B, B > E : h(s) = n0 6 ^ 9 I — - E 2 VB 2 -E 2 cot es^B 2 - E 2 ] . 

B 2 — E 2 \es J 

(13) 

In the cases I and II h(s) has no zeros for real s so one kan keep the contour as 
in Eq.(H). A special case of III is when E — > and we get back the result in M by 
deforming the contour to the imaginary axis. Another interesting special case of E _L B 
is when E = B in which case we have a = b = 0, h(s) = 1/s, just like in absence of 
the external field, but with a potential A . We then find 

C^(F, U = 0, Ao) = C^(E ±B,E = B) 

dp f F (p , A Q )6((p - eA ) 2 - m 2 )((p - eA ) 2 - m 2 f 2 , (14) 



3tt 2 

a result similar to the absence of quantum corrections in a propagating plane wave 
at zero temperature 0. Notice that, because of the lack of Lorentz invariance, this 
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observation is non-trivial since C^r does not only depend on T and Q. To obtain 
Eq.fln|) we have used dimensional regularization as in |7[ and the relation 



oo n 

ds s v sinfas J 

o 4 7 



if a > 



sinf(z/ + |) ifo<0. 



sm + ^ 

3 
2 



(15) 



3 Weak E— field expansion 



The ^-function in Eq.(|T4|) arises from the non-analytic behaviour of Eq. (|15"D . If we 
take B = and expand the integrand of Eq.(|7|) in powers of E, using Eg-flHf), we 
find that the s-integral can be performed term by term using Eq.(^), resulting in the 
same ^-functions. However, the remaining p -mtegral becomes infra-red divergent at 
(Po ~ ey 4o) 2 — m2 - This problem occurs already at 0(E 2 ) in contrast to the case of 
a B-field 0] where the 0(B 2 ) term can be calculated by a direct expansion of the 
integrand. It is not a priori clear that a power expansion in E 2 exists, and it is also 
plausible that the expansion is only asymptotic just like at zero temperature. Such 
an expansion would still be useful in many physical situations where the field is not 
too strong. To find the leading weak field expansion we have to be more careful than 
simply expanding the integrand. Using the notation from Eqs.(|H||n]), we derive from 



£&W-£5?{0) = -^r oD diufF{uJ + eAo) Io 



ds 



(h(s)s) 1/2 cos \uo 2 s 



l-h(s)s 
h(s)s 



7T 



cos ( {u — m )s 



(1/2 :.. / .2 1 h(s)S \ I , 2 ; :! , 



-(h(s)sY /2 sin [uj 2 s- 



71 



sm (a; — m )s — 



h(s)s J ~~ V" "" '' A, t 
= C l (E)+C 2 (E)] ' (16) 

where u = po — eA . In the first part of the curly bracket in Eq. flT6|) , called Ci, we can 
use the expansion h(s)s ~ 1 + (seE) 2 /3 to find the finite E 2 term 

{eEf r 00 duo 



C X {E) 



—f F (u + eA )9(uj 2 - 



m 



(17) 



24tt 2 7-oo (oo 2 -m 2 ) 1 / 2 
The infra-red problems arise when trying to expand £2 in E. After introducing a new 
integration variable x = u 2 , and performing a partial integration with respect to x, we 
get 

1 



= g^y Q dxJ Q -^(h(s)s)V 2 cos(xs-m 2 s-^) 



d 

dx 



f F (u + eA ) + f F (-u + eA ) . / 1 - h(s)s" 
sm xs- 



h(s)s 



(18) 
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which can be expanded to 0(E 2 ). To this order in the field we then have that 



^effl ) ^effW 247r2 J m ( u 2_ m 2y/2 du \ e p(« J +A -ri e ^-A +») + 1 

Eq.(|l9D leads to 



(19) 



= ^X(O) - V ' /(1/i-eAol-m) , (20) 

V (A* — e A)) — 



if T = and 

= #0) - ^ (l - <)- g ( (m/r)2 + ^ + l)vf2 ) . (2D 

if — eA = 0. The Debye mass can be extracted as the second derivative of C^{E) 
with respect to A . From the zero field part in Eq.([TJp we get m 2 (/x) ~ e 2 (/i — eA ) 2 /ir 2 
and m 2 (T) ~ e 2 T 2 /3 in the high density and temperature limit. This agrees with 
[ |T3| , [TJj] but not with |[L1|| . Corrections from finite .E-field, to lowest order in the field, 
can also be calculated making use of Eq. (p~9|) and we find for high temperature and 
high density, respectively 



., , a (eE) 2 31 t/ , m 4 / „„ / ' m,., , , 

m 2 (T,E) = mi( r)--^^C(5)j 7 (l + 0((^]) • (22 



where £(n) is the Riemann zeta-function and a = e 2 /47r is the fine-structure constant. 
We notice that the E'-field tends to decrease the screening mass. 

In analogy with WB, [/J we can define an effective fine structure constant by 

111 dC^(E) 



(23) 



a(T,fj) a aE dE 

In the high density or temperature limit we then have that a(T,fi) — > cx/(l — a/37r) 
showing a completely different behaviour from the a defined using a magnetic field pj 
which satisfies a zero temperature renormalization group equation. 



4 Pair production 



One issue that has been discussed in the literature is pair production at finite T [[H], [11 
which we find to be absent from the thermal part of the one-loop contribution, in the 
real-time formalism. This can be seen directly from Eq.(|l]) in which the thermal part 
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is manifestly real. It is illustrative to compare with the standard calculation of the 
effective potential V(4>) for a spontaneously broken A0 4 -model where one may find 
an imaginary part when the effective mass, M 2 (0) = A0 2 /2 — m 2 , is negative. The 
expression 

dV(<f>) _\<p f d 4 p ( % % \ 

# ~ 2 J (27r) 4/ ^ oj \p 2 — M 2 + ie p 2 - M 2 - te ) ' ^ ' 

is obviously real for any real M 2 . An imaginary part of V(<f>) is found when first calcu- 
lating Eq. ( ^4]) for a positive M 2 in the limit e — > and then performing an (ambiguous) 
analytic continuation in M 2 . This new function corresponds to, for negative M 2 , a po~ 
integration contour which does not follow the real axis but goes above (below) the poles 
for negative (positive) p , or vice versa. The procedure of obtaining an imaginary part 
in this way is thus both ambiguous and does not correspond to the e-prescription in 
Eq.fl24p. The conclusion we can draw from the example above is that the standard 
real-time calculation of the one-loop effective potential gives no imaginary part. One 
should, however, remember that the standard real-time rules are usually derived under 



the assumptions of certain factorization properties (see JT6| for a discussion) that may 
not be fulfilled in presence of unstable modes. An imaginary part can be obtained in 
the imaginary-time formalism, but it is not clear that a consistent calculation of the 
pair production rate can be performed assuming equilibrium. 

We conclude with a remark on heavy quark production in ultra-relativistic nuclear 
collisions, showing the physical importance of pair production. This example serves as 
a motivation for finding a solution to the problem with the imaginary part. If JZ g -> q q 
is the production rate of thermal gluon decay into quark-antiquarks one finds e.g. that 



17 



%?-0.01 , (25) 

for quark masses m such that m/T < 1. Taking the imaginary part of the effective 
action ImL^fJE) found in [IT, [J seriously, which we question in our paper, we would 
find for T>m and /i = 0a rate 



ImC'fiW ( 0.3 



T 4 \T, 



- GeV 



(26) 



which can be as important as perturbative production rates. Here we have used that 
eE ~ lGeV 2 . It turns out that perturbative production of heavy quark-antiquarks 
due to the time-variation in the background electric field can be as large as the ther- 
mal gluon decay rate Eq.(p5|). Since quark-antiquarks are spontaneously produced by 
Schwinger's mechanism, the background electric field will be depleted, as mentioned in 
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the introduction. The time-scale, ta, for the depletion of the electric field we estimate 
using Schwinger's expression for the rate of pair production, w s , and mtdW s ~ E 2 /2. 
Approximatively, this amounts to assuming an exponential decay of the background 
electric field, i.e. eE(t) ~ eEexp(— t/tA. The explicit time-dependence in the elec- 
tric field leads to pair production already at the one-loop level with a rate TZe- A 
straightforward calculation leads to the result 



where we have used a s = 0.32 and considered a quark mass m = 0.2 GeV. If the 
pair-produced particles are assumed to be in thermal equilibrium, one can expect a 
decrease of this production rate from Pauli blocking factors. 
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